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Abstract: We consider jet-shape observables of the type proposed recently |1T|, H 
where the shapes of one or more high-p-r jets, produced in a multi-jet event with defi- 
nite jet multiplicity, may be measured leaving other jets in the event unmeasured. We 
point out the structure of the full next-to-leading logarithmic resummation specif- 
ically including resummation of non-global logarithms in the leading- iVc limit and 
emphasising their properties. We also point out differences between jet algorithms 
in the context of soft gluon resummation for such observables. 
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1. Introduction 

It has long been known that probing the shape and structure of high-p^ jets is poten- 
tially of great value in searches for new particles at collider experiments . With the 
advent of the LHC and much activity in improving and developing jet algorithms 0- 
studies of this nature have received considerable impetus. In particular, much 
recent attention has been focused on using jet studies for the identification of boosted 
massive particles which decay to hadrons forming a coUimated jet, see for instance 
Refs m-l^- 

In the same context a method has been recently suggested to study the shapes 
of one or more jets produced in multi-jet events at fixed jet multiplicity |I[, 0. The 
precise details of the observable suggested in those references involve defining a jet- 
shape energy- fiow correlation similar to that introduced in Ref. [0. Specifically the 



proposal was to measure the shapes of one or more jets in an event leaving other jets 
unmeasured and introducing a cut on hadronic activity outside high-pT (hard) jets, 
to hold the hard-jet multiplicity fixed. This is in contrast to for instance hadronic 
event shapes ^ ^ which by construction are sensitive to the shape of the 



overall event rather than an isolated jet. 
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In the present paper we wish to use this observable study to make 

several points that we believe will be useful both within and outside the specific 
context. The main points that we wish to address concern the resummation of soft 
gluon effects that become important in describing the observable distribution for 
small values of the shape variable p and the energy cut Eq. 

In particular in this paper we shall address the structure of non-global logarithms 

^ as well as compute them in the large- iVc limit for jets defined in the anti-kt 
algorithm. We remind the reader that observables that are sensitive to radiation in 
a limited phase-space region such as the interior of a particular jet are non-global in 
the sense that they receive logarithmic contributions from correlated soft emission, 
which are highly non-trivial to treat to all-orders. Existing resummations of non- 
global logarithms have been confined to a few special cases |3^, ^ ^7\, ^ and to 
the large- A'^c limit. Given that the observable we study in the current paper is non- 
global, it is worth examining in detail the precise structure of non-global logarithms, 
which by definition start at Ola^) in the soft function and are of the same size as 
the logarithms resummed in Refs. and thus need to be considered as well in 
order to achieve NLL accuracy. In the current paper we find non-global logarithms 
arise both in the ratio of the the energy cut-off Eq and the shape variable p as well 
as in Eq/Q where Q is the hard scale of the process, which is naturally of the order 
of the hard jet pt^. More to the point we argue that in the limit of narrow well- 
separated jets a simple picture emerges for non-global logarithms. The simplicity 
in the non-global structure is to do with the fact of QCD coherence. Narrow well 
separated jets do not affect each others evolution even in the non-global component 
which arises individually as an edge effect from the boundary of each jet, precisely 
as the non-global logs in the case of a hemisphere mass in e^e~ annihilation arise 
from the edge separating the observed and unobserved hemispheres [^. Hence the 
resummation of non-global logarithms arising at each jet boundary can simply be 
taken from the existing result for a hemisphere^ up to corrections that vanish as 
powers of the jet radius. The simple structure of non-global effects in turn provides 
us with an ansatz that can be used for any jet event of arbitrary jet multiplicity. 

We also assess here the numerical contribution of the non-global logarithms and 
find that while limiting the value of Eq is of some use in diminishing their size the 
effect is still of order twenty percent as far as the peak height of shape distributions is 
concerned. In fact we find that changing the value of Eq is not particularly useful as 
a means of reducing the non-global contribution. Specifically following the original 
proposal in Ref. |^ it was suggested in Ref. that one may take the value of Eq/Q 
to be of the same order as the jet shape variable p, which we agree eliminates the 
non-global contribution from the measured jet. However in this case the contribution 



^The potential presence of such logarithms was also mentioned in [Q, ^. 
^This statement should be qualified as it is correct only for the case of the anti-fcj jet algorithm | 
which is the one we recommend for study of such observables. 
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from the unmeasured jet becomes as significant as the contribution we are attempting 
to ehminate and hence the overall effect of this choice turns out to actually increase 
the overall non-global component. With the resummed results of the current article 
however one does not have to be too concerned about the precise choice of Eq as the 
non-global terms should be accounted for, at least within the large- iVc approximation 
and up to corrections vanishing as powers of the jet radius ^. 

Another point we wish to make is concerning the role of the jet algorithms. 
The computation of non-global logs in the leading- iVc limit can actually be carried 
out in any jet algorithm by means of the numerical codes developed for instance in 

p9t 0. Indeed it was found |^ that the use of certain sequential recombination 
algorithms (such as the kf or Cambridge- Aachen (C-A)) can significantly reduce the 
non-global logarithms due to the soft gluon clustering inherent in such algorithms. It 
was however later demonstrated that one pays a price for this reduction in the 

non-global component in the form of a more complicated result for the independent 
emission terms. While independent emission is commonly associated with the expo- 
nentiation of the single-gluon result, this association is spoiled by the application of 
sequential recombination algorithms other than the anti-/cj algorithm. As we shall 
show in this paper the result of soft gluon clustering in the kt and C-A algorithms 
modifies the independent emission (global term) which deviates from the naive ex- 
ponentiation of a single gluon at a relevant single logarithmic accuracy. Moreover 
the effect of the clustering near the boundary of a collinear jet no longer produces 
logarithms suppressed in the jet radius R as was the case for small central rapidity 
gaps discussed in Refs. but rather pure single-logarithms independent of R. 

These effects are absent for the anti-kt algorithm as already pointed out in Ref. [^, 
since that algorithm clusters soft gluons independently to the hard parton and hence 
produces circular jets in the soft limit, i.e. it can be regarded in this limit as a rigid 
cone. Hence for the present moment and pending a resummation of the clustering 
logarithms along the lines of that carried out for gaps between jets |Q we confine 
our studies to the anti-fct algorithm. We do however provide an explicit fixed-order 
computation of the single-logarithmic corrections in the independent emission piece, 
that arise in other algorithms as we believe this point deserves some stress. 

The paper is organised as follows. In Section 2 we define our observable, choosing 
the jet mass in dijet events as an example of a jet-shape observable, while imposing 
a cut Eq in the inter-jet energy flow. In Section 3 we perform the leading and next- 
to leading order calculation of such observable in the soft limit, which elucidates 
the structure of the logarithms arising from independent soft gluon emissions as 
well as non-global logarithms from correlated emissions. We use these results to 
construct an argument which culminates with the resummation of these logarithms 
in Section 4. We also present a study which assesses the numerical significance of 

■'These effects would amount to perhaps a ten percent change in the non-global term which we 
do not expect to be of significant phenomenological consequence. 
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the non-global logarithms as a function of the parameters p and Eq. In Section 5 we 
highlight the fact that for algorithms other than the anti-fct exponentiation of the 
single gluon result is not sufficient to capture the next-to-leading logarithms even in 
the independent emission piece, let alone the non-global terms. Finally we draw our 
conclusions in Section 6. 



2. High-pT jet shapes and inter-jet energy flow 

We wish to examine a situation where one studies the shapes of one or more high-p-p 
jets in jet events with definite multiplicity. From the results we shall obtain below 
for such events one can draw conclusions also about the single inclusive jet mass 
distribution for instance for the process pp ^ j + X, where one can demand the 
production of a jet j setting a value for a particular jet-shape, while summing over 
everything else denoted by X. 

For the points we wish to make in this paper we can for illustrative purposes and 
without loss of generality consider high-p-r dijet events. In order to restrict the jet 
multiplicity we can place a cut Eq whereby we veto the inter-jet activity such that 
the sum of transverse energies of emissions outside the two high-pT jets is less than 
this value. This was also the definition adopted in Refs. [jl|, |^ where the parameter 
A indicated a cut on additional jet activity along the above lines. 

Moreover, in this paper we are interested in physics at the boundary of the 
triggered hard jets and specifically in the non-global logarithms that arise at these 
boundaries. Hence we can for our discussion ignore the effects of initial state radiation 
which can simply be accommodated once the structure of the results is understood. 
Since it is this structure we wish to focus on, it proves advantageous to consider as an 
analogy the production of dijets in e+e~ annihilation which enables us to ignore the 
detail of initial state radiation. Hence all our points can be made in full generality 
by considering two hard jets in e~^e^ processes where one measures the shape of one 
of the jets leaving the other jet unmeasured as prescribed in Refs. [|I], 0]. Our results 
should also then be directly comparable to those obtained by other authors using 
soft-coUinear effective theory 



2.1 Observable definition 

We shall pick the jet mass as a specific simple example of a jet-shape variable though 
one can consider also, for instance, the angularities first studied in The 
observable we study has the same logarithmic structure as the distribution in the 
angularity corresponding to a = 0. We study the shape cross-section 
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where ao is the Born cross-section and p'l denotes the normahsed jet-mass of the jet 
with momentum Pi. The above equation indicates that we are restricting the mass 
of the jet with three momentum Pi to be less than p leaving the shape of the other 
jet with momentum P2 unmeasured. We have also restricted the inter-jet energy 
flow Eq to be less than Eq as discussed and hence our observable definition above is 
in precise accordance with the definition in Refs. ||T], 0. We shall in future leave the 
dependence on jet momenta Pi, P2 unspecified and to be understood. 

We wish to carry out a calculation for the above observable which includes a 
resummation of large logarithms in R"^/ p to next-to-leading logarithmic (equiva- 
lently single logarithmic) accuracy in the exponent. We include a description of 
non-global single logarithms in the leading- A^^c limit. Additionally we wish to resum 
the logarithmic dependence on Q/i?o, where Q is the hard scale of the process, to 
single logarithmic accuracy again accounting for the non-global contributions. Our 
main aim is to study the effect of the non-global logarithms neglected for instance 
in previous calculations of jet shapes [|T], ^ on the cross-section Eq. (^.11). While 
resumming logarithms in p and Eq we shall neglect those logarithms that are sup- 
pressed by powers of the jet radius R which shall enable us to treat non-global 
logarithms straightforwardly Hence our calculation addresses the range of study 
where Eq/Q ^ p and is valid in the limit of relatively small jet radius R. We shall 
not however resum terms varying purely as ag In R which for the values of R we 
consider can safely be ignored from a phenomeno logical viewpoint. We thus aim to 
resum large logarithms in p and Eq/Q in what one may call the approximation of 
narrow well separated jets. According to our estimates this approximation and our 
consequent resummation should enable relatively accurate phenomenological studies 
of jet shapes. 

We shall begin by carrying out a calculation of the logarithmic structure that 
emerges at the one and two gluon levels, in the limit of soft gluon emission. These 
calculations will help us to identify the full logarithmic structure and point the way 
towards a resummed treatment. We start below with a leading order calculation in 
the soft limit. 



3. Soft limit calculations 



We start by considering the effect of a single soft emission by a hard qq pair, produced 
in e+e~ annihilation. At this level all infrared and coUinear (IRC) safe jet algorithms 

^More specifically we shall neglect corrections varying as R'^/Aij where Aij = 1 — cos 9ij is a 
measure of the angular separation between the hard jets. This parameter emerges naturally in 
fixed-order computation of non-global logarithms for energy flow outside jets |Q and it was also 
treated as negligible in [0, ||. 
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will yield the same result. We can write the parton momenta as 



Pi 



0,0,1) 
§(1,0,0,-1) 



(3.2) 
(3.3) 



(3.1) 



k 



u {1, sin 9 cos (f), sin 9 sin (p, cos 9) 



where pi and p2 are the hard partons and we have neglected recoil against the soft 
gluon emission k, which is irrelevant at the logarithmic accuracy we seek. Let us take 
the momentum pi to correspond to the measured jet direction. Hence if the gluon is 
combined with the parton pi one restricts the mass of the resulting jet to be below 
p while if combined with p2 the mass is unrestricted. Likewise one can consider the 
parton p2 to be in the measured jet direction, which will give an identical result. 

Introducing the jet mass variable p — Mj/Ej, where one normahses the squared 
invariant mass M| to the jet energy Ej one can write 



where in our soft approximation the jet energy is set at Q/2. Note that there is no 
constrain on the gluon energy when it is combined with the jet j2. 

In all the commonly used IRC safe jet algorithms the soft gluon k will form a 
jet with a hard parton if it is within a specified distance R of the hard parton. The 
distance is measured for hadron collider processes in the (77, 0) plane as Ar]'^ + 
where At] is the separation in rapidity and A0 is the separation in between the 
hard parton and the gluon k. In the limit of small angles, relevant for small R values 
R <^1, which we shall consider here, the distance measure reduces to ^^j. where 9pk 
is the angle between the gluon k and hard parton p. Thus k and pi form a jet if 
^lik < Otherwise the gluon is outside the jet formed by pi which at this order 
remains massless. If the gluon does not also combine with hard parton p2 to form 
a jet, one restricts its energy to be less than Eq as required by the definition of the 
observable. Differences between the various algorithms shall emerge in the following 
section where we examine the emission of two soft gluons. 

Thus we can write for the contribution of the real soft gluon k with 9^^,^ < R^ 



where we restricted the jet-mass to be less than the specified value pQ^/4: and the 
superscript r denotes the real emission piece. In this same soft region virtual correc- 
tions are exactly minus the real contributions, but unconstrained; therefore we can 
cancel the real emission result above entirely against the virtual piece and we are left 




(3.4) 




(3.5) 



6- 



with 



/ / R2^[-Fr-p] ^ (3.6) 

where we constructed the jet mass p = 4^(1 — cos^^) ~ and where we used 

the small- angle approximation since 6^ < <^ 1. The suffix "in" denotes the 
contribution to S from the region where the gluon is in the measured jet. Performing 
the integral over angle with the specified constraint results in 

E,. = -£i^ r U^-f] = 5!e [ff - p) . (3.7) 

JpQ/2R^ ^ \ QP J 27r p ^ 

Next we consider the region where the soft emission flies outside either hard 
jets, with the corresponding contribution Sout- Since here we are no longer confined 
to the small angle approximation we use kt and 7] with respect to the jet axis as 
integration variables where r] is the gluon rapidity. In these terms one can represent 
the contribution of the gluon k after real-virtual cancellation as 



2CFas f dkt /-i^^/fl 



-'out 



TT J kt J_ 



dr]Q {kt cosh T] — Eq) . (3. 

In2/i? 



where the limits on the rapidity integral reflect the out of jet region. Performing the 
integrals we get to the required single-logarithmic accuracy 

So,t = -2Cf^ In (Q/Eo) (^2 In |^ . (3.9) 

The full soft result at leading order is Si = + Sout- As is well known the jet- 
mass distribution receives double logarithmic corrections which in the present case 
are in the ratio R^/p- Taking account of hard coUinear emissions one would obtain 
also single logarithms in R^/p, which we shall account for in our final results. 

The above calculation having set the scene we shall now move to considering 
two-gluon emission and the structure of the non-global logarithms that arise at this 
level. 

3.1 Two-gluon calculation and non-global logarithms 

Going beyond a single soft emission to the two gluon emission case the precise de- 
tails of the jet algorithm start to become important. In what follows below we shall 
consider only the anti-fcj algorithm since in the soft limit the algorithm functions 
essentially as a perfect cone algorithm ^ . In particular this implies that soft gluons 
are recombined with the hard partons independently of one another (one can neglect 
soft gluon clustering effects) which considerably eases the path to a resummed pre- 
diction. The logarithmic structure for other jet algorithms is also interesting and we 
shall discuss it in a later section. 
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Here we carry out an explicit two-gluon calculation to obtain the structure of non- 
global logarithms for the observable at hand. Referring to the non-global contribution 
to E as S, we compute below 5*2 the first non-trivial term of S. Our results shall 
indicate a way forward towards a resummed result incorporating these effects. As in 



Refs. ^ we shall consider the emission of gluons ki and k2 such that coi ^ 0J2-, 
i.e. strong energy ordering. In this limit the squared matrix element can be split into 
an independent emission term oc C"^ and a correlated emission term oc CfCa The 
former is incorporated in the standard resummed results based on exponentiation of 
a single gluon, which we discuss later. 

Let us concentrate on the CfCa term missed by the single gluon exponentiation, 
and which generates the non-global logarithms we wish to study and resum. We now 
consider the following kinematics: 

pi = |(l,0,0,l) (3.10) 

P2 = |(l,0,0,-1) 

ki = ui (1, sin 6*1, 0, cos6'i) 

^2 = 1^2 (1, sin 6*2 cos 0, sin 6*2 sin 0, cos 6*2) 

(3.11) 

with ui ^ U2- 

Let us consider the situation where the harder gluon ki is not recombined with 
either jet but the softest emission /c2 is recombined with pi. This situation corre- 
sponds to the diagram of the left in Fig. |I[ In the small-angle limit, which applies 
for the case i? <^ 1, the condition for ^2 to be recombined with pi is simply 6*1 < 
or equivalently 1 - < cos 612 < 1 while one has -1 + < cos < 1 - R^/2 
which ensures that ki is outside the jets. We integrate the squared matrix element for 



ordered soft emission |^ over the azimuth of gluon ^2 to get the angular function p5 



^ = ("3 12) 

(cos6'2 - cos6'i) (1 - cos6'i) (1 -Fcos6'2)' 

Then defining the energy fractions Xj = the required integral for the non-global 
logs reads 

2 dxi dx2 „ f 2Eq . ^ . . 

- xi \<S(xi- X2) 



S2 = -4CFC4^y [ dj±f dX2Q 

\2nJ Jq Xi Jq X2 



Q 

I d COS 62 I d cos ^1^9 (2x2(1- cos ^2) -p), (3.13) 

il-_R2/2 J-l+_R2/2 



where we note the constraints on ki and ^2 imposed by the observable definition. 
Note that as in Ref. [B^ the constraint on ^2 emerges after including the term where 
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^2 is a virtual gluon such that the divergence of real emission is cancelled and the 
piece we retain above is the virtual leftover. Integrating over xi,X2 we obtain 



/a \2 /-i rl-RV2 
S2 = -2CfCa(^) / rfcos^a / ^cos^i 



2 pQ (, pQ 



"' 4E„(r-co..,) n'" 4^.(l-cos.J "- f^"' 

The angular integrations over 6*1, 6*2 gives the number 7r^/6 provided we neglect 
terms of order E? and those varying as pQ/{2EqR^). We recall that as stated before 
we neglect logarithms suppressed by powers of R and also that our resummation will 
be valid when p/ <^ Eq/Q and hence can ignore the corrections to 7r^/6. 

Thus in the small R limit the result for the leading non-global piece is 

S2 = -CfCa[-) -In -^^[-^-P)- (3-15) 

We note that one also can receive a contribution to the non-global logs from 
the case where ki is part of the unmeasured jet but this configuration produces a 
coefficient that varies as R^ and hence can be ignored, consistently with our approx- 
imation. Lastly carrying out the integration with the harder gluon ki inside the 
measured jet and the softest one ^2 outside does not give us large logarithms in the 
region we are interested in, hence Eq. (|3.15| ) is our final result for the first non-global 
piece affecting the p distribution. 

Next we consider the case that the harder gluon ki is in the unobserved jet and 
emits k2 outside both jets, as depicted in Fig. |1|, on the right. In this case repeating 
the calculation in the same way produces to our accuracy 

S2 = -C,cA^)"^-\n'^. (3.16) 
V27r/ 3 2^0 ^ ' 



The above results are noteworthy in many respects. Note that the result Eq. ( p.l5| ) 



corresponds to the result already obtained for the hemisphere jet mass in Ref. p5| 
provided one replaces 1/p in that result by ^^/-^- This is because the non-global 
evolution takes place from energies of order Qp/R^ up to those of order Eq, whereas 
for the hemisphere mass the relevant energy for the harder gluon was of order Q. 
More interestingly the coefficient of 5*2, vr^/3, is the same as was obtained there. The 
origin of this is the fact that the coUinear singularity between ki and ^2 dominates the 
angular integral. As has been noted before |Q as one separates the gluons in rapidity 



the contribution to the non-global term, which represents correlated gluon emission, 
falls exponentially as gluons widely separated in rapidity are emitted essentially in- 
dependently. Thus in our present case, up to corrections suppressed by R^ the results 
for the p distribution arise from the edge of the measured jet independently of the 
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Figure 1: Diagrams representing the correlated emissions which give rise to the lowest- 
order non-global logarithms. On the left: the harder gluon ki lies outside both jets and 
the softest one k2 is recombined with the measured jet and contributes to the jet-mass 
distribution. On the right: the harder gluon is inside the unmeasured jet and emits a 
softer gluon outside both jets, which contributes to the iiJo-distribution. 



evolution of the unobserved jet. Likewise there are non-global logarithms given by 
Eq. ( |3.16| ) which affect purely the inter-jet energy flow Eq distribution. These arise 
purely from the edge of the unmeasured jets and are independent of the evolution of 
the measured jet which is well separated in rapidity (similar results were obtained in 
the work of Refs. [p| p|l). 

Thus a simple picture arises for non-global logarithms where each jet evolves 
independently and the effects arise from the edges of the jet with logarithms involving 
the ratio of the shape variable p/R"^ to the energy flow variable Eq/Q coming from 
measured jets and unmeasured jets independently contributing logarithms in Q /Eq. 
The coefficients of the non-global logarithms will be identical within our accuracy 
to those computed for the hemisphere mass (where the effect is again an edge effect 
coming from the hemisphere boundary) and hence the resummation of the non- 
global effects from each jet can simply be taken from the resummation carried out 
in Ref. simply modifying the evolution variable. This will be done in the next 
section. 

To conclude we wish to draw attention to the fact that we have determined, with 
a fixed order calculation, the precise non-global structure which was not included in 
Refs. 10, H and knowledge of the nature of these logarithms should pave the way 
for more accurate phenomenological studies. We remind the reader that our study 
above is valid only for the case of the anti-Zcj algorithm. Other jet algorithms will 
give different non-global pieces as discussed in Refs. |^0|, ^ In fact even the 
resummation of independent emission terms will be different in other algorithms, a 
fact that is not widely appreciated and that we shall stress in a later section. 

In the following section we turn to resummed results and provide a simple ansatz 
which will be valid for arbitrarily complex processes involving jet production. 
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4. Resummation 



Having observed the key feature of the non-global logarithms (independent contri- 
butions from each jet) that allow us to write a resummed result we shall now focus 
on the resummation in more detail. The main point to note is that the non-global 
logarithms provide a factor that corrects straightforward single-gluon exponentia- 
tion 

j:(^Q]= ( ^, R-] s-^ f # V (4.1) 



p ' Eq) \ p ' EqJ \ Qp ' Eo^ 

Thus we shall first provide the result for the single-gluon exponentiation taking 
account of hard-coUinear emission and the running coupling, which contains leading 
and next-to-leading logarithms in R'^/p as well as leading logarithms in Q/Eq. 

4.1 Independent emission contribution 

The resummation of independent emission contributions based on a squared matrix 
element that has a factorized structure for multi-gluon emission is by now a standard 
procedure and we shall avoid listing these details (see for instance for a detailed 



study of these techniques). We shall provide here only details of the final result 
for independent emission valid for the anti-fct algorithm only. We stress once again 
that even the independent emission piece will differ at next-to-leading logarithmic 
accuracy from that reported below if using another jet algorithm. 

The result for the independent emission contribution can be written in the usual 



form 47 



\p ^(l + 7^g """^ ^ ' 

Here TZp and TZeq are functions of I p and Q / Eq respectively, representing the 
exponentiation of the one gluon result. They describe the resummation of large 
logarithms to next-to-leading logarithmic accuracy in R^ / p and leading logarithmic 
accuracy in Q/Eq except for the inclusion of non-global logarithms not described by 
independent emission of soft gluons. 

With inclusion of running coupling effects and the effects of hard coUinear emis- 
sion the function TZp can be written as 



= V / ^^^^sihmki) , (4.3) 



where we defined 



3 



+ H■^)e(^llefI-i^,V ,4.4) 



pQ J \A q^) 4/?2 
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where the factor e~'^/^ in the argument of the logarithm in the first term above takes 
account of the hard colhnear region 2a;/Q — )■ 1 ^. 

Carrying out the integral over kt one obtains the familiar result for TZp as follows 

7^, = -L/l(A)-/2(A), (4.5) 

and 



d_ 
di 

The functions /i and /2 are listed below 



^; = -— • (4-6) 



/i(A) = [(1 - 2A) In (1 - 2A) - 2 (1 - A) In (1 - A)] , (4.7) 



and 



/.(A)^-^[21n(l-A)-ln(l-2A)]-^ln(l-A) 



In (1 - 2A) - 2 In (1 - A) + ^ In^ (1 - 2A) - In^ (1 - A) 



(4i 



A = /SoagL, L = In and = as {QR/2) is the MS strong coupling. In the above 



results the /3 function coefficients /3o and /3i are defined as 

IICa - 2nj nC\ - hCAUf - 3C^n^ 
12^^ ' 24^^^ ' ^^-^^ 



and the constant K is given by [48 



/67 



TT 



Likewise for the function TZeq we have 

2C 2 

nE, = ^ln-ln(l-2A), (4.11) 

where here A = PoOgL, L = In and as = as (<5/2). Note that the function TZp 
contains both a leading logarithmic term L/i(A) and a next-to-leading or single log- 
arithmic term /2(A) while the leading logarithms in TZeq are single logarithms and 
next-to-leading logarithms in this piece are beyond our control. The term F (l + TZ'p) 
arises as a result of multiple emissions contributing to a given value of the jet-mass 
and is purely single-logarithmic. The corresponding function for the Eq resumma- 
tion would be beyond our accuracy and hence is not included. We note the results 
presented here for TZp are identical to the ones for the e^e~ hemisphere jet-mass [^7]], 
with the replacement p — > p/R^ and as{Q) — > as{QR/2). 

^In order to obtain this one replaces as usual ^ dx where x — 2(jj/Q in the integral 

over gluon energy, which is essentially introducing the full splitting function instead of its soft 
singular term. 
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4.2 Non-global component 

The non-global terms arise independently from the boundary of individual jets in 
the approximation of narrow well-separated jets. The result for an individual jet 
is the same as that for energy flow into a semi-infinite rapidity interval which was 
computed in the large- A'^c limit in Ref. 

In our two-jet example the contribution of non-global logarithms can be thus be 
written as 

measured) "S* (tujimeasured) ; (^■-^^) 



where the function S{t) was computed in Ref. p5|. From that reference one notes 
that 

SW^exp(-C.C.4(i±M!)*^), (4.13) 

where a = 0.85C^, b = OMCa, c = 1.33. 

The single logarithmic evolution variables for the measured and unmeasured jet 
contributions read 

J- dx 

Wasured = 7^ / „ as{xEo) , (4.14) 

cc 

52- 









/ pQ 




2EqR 












Q 



^ dx 

^unmeasured / Ois {•^Q/'^) j (^■-^^) 

2tt J'iEiL X 

Q 

which represent the evolution of the softest gluon with a running coupling that 
depends on the gluon energy. For the measured jet the softest gluon evolves between 
scales of order Qp/ and Eq while for the unmeasured jet the evolution is from Eq 
up to the jet energy Q/2. 

Carrying out the integrals (one-loop running coupling is sufficient here) gives 

1 / 2E R? \ 

^measured = -^—T In I 1 - /^otts (-^o) In — ^ I , (4.16) 

4vr/?o V QP J 

^unmeasured = hi ( 1 - /^otts (<5/2) lu ) . (4.17) 

In the following sub-section we shall illustrate the effects of the non-global logarithms 
on the shape- variable distributions for different values of Eq. 

4.3 Numerical studies 

Let us examine the impact of non-global logarithms on the differential jet mass 
distribution, divided by the inclusive rate; at our level of accuracy we have: 

1 do" dS 

-— = — , 4.18 
cr dp dp 

with S given by Eq. ( ^4.1| ). From Figure. ^ one can see that non-global logarithms 
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Figure 2: The jet mass distribution ~ for Q = 500 GeV, R = 0.4 and Eq = 15GeV 
(left) and 60 GeV (right). The curve in dotted red corresponds to neglecting non-global 
effects while that in solid blue takes them into account. 



do not change significantly the position of the peak of the distribution. However, 
their inclusion leads to a reduction in the peak height of fifteen percent or so for 
Eq = 15 GeV and about twenty percent or so for Eq = 60 GeV. Increasing Eq 
further one will observe that the effect of non-global logarithms on the peak height 
can be as significant as about 30%. The plots above are for Q = 500 GeV which may 
be translated into a jet pt of about 250 GeV or so at a hadron collider. 

On the other hand it has been suggested that one may eliminate non-global 
logs by choosing Eq/Q of order p. In our case (small R) this prescription amounts 
to the choice 2EqR^ = Qp. While this rids us of non-global logarithms from the 
observed jet boundary the contribution from the unobserved jet boundary becomes 
increasingly important. This is reflected in Figure |^ which plots separately the 
factors 1 — 5'(tmeasured) (dotted red curve), 1 — 5'(tunmeasurcd) (dashed blue curve) and 
1 — S"'^ (solid green curve), where S"^ is defined in Eq. ( [4.12|) as the product of the 
S factors . The plots are presented as functions of Eq for the illustrative value of 
p = 5 X 10~^. Other parameters are the same as for the previous plots. As one 
can readily observe increasing the value of Eq leads to a growth of the non-global 
contribution from the measured jet while the contribution from the unmeasured jet is 
somewhat diminished. Lowering Eq leads to the opposite effect and the unmeasured 
jet contributions become increasingly significant. It is noteworthy that changing the 
value of Eq in the range indicated has no significant effect on the size of the non- 
global effect overall. Also worth noting however is that the choice £"0 = pQ/{2R^) 
(the lowest value of Eq shown in the above mentioned plot) which eliminates the 
contribution from the measured jet (i.e. the red curve goes to zero) is not very 
helpful as the overall contribution stemming from the unmeasured jet entirely is 
more significant than for the higher values of Eq discussed before. From this one 
realises that progressively decreasing the value of Eq is not a way to eliminate the 
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Figure 3: Non global contribution 1 — S{t) from the measured jet (dotted red), the 
unmeasured jet (dashed blue) and overall (solid green) as a function of Eq for p = 5 x 
10-^, R = 0.4. 



non-global contribution, for the observable at hand. 

In the following section we shall show that for algorithms other than the anti-fct 
even the independent emission resummed result is not equivalent at next-to-leading 
logarithmic level to the exponentiation of the single-gluon result. 



5. Other jet algorithms 

Let us now consider the situation in other jet algorithms where the clustering or 
recombination of soft gluons amongst themselves may be an important effect. One 
such algorithm is the inclusive kt algorithm discussed for the case of central gaps 
between jets in Refs. For such algorithms, starting from the two-gluon 

level, we need to revisit the independent emission calculations and correct the naive 



exponentiation of a single gluon. Note that in Refs. [^2|, ^ the single logarithms 
obtained as a result of clustering were proportional to powers of the jet radius which 
would make them beyond our control here. However, as we shall see, in the collinear 
region we are concerned with here, this power suppression does not emerge, making 
these logarithms relevant to our study. To illustrate the role of soft gluon clustering 
and recombination we focus on the on the single inclusive jet-mass distribution and 
we ignore the cut corresponding to Eq. Placing this cut does not affect the conclusions 
we draw here. 

To set the scene let us first carry out the independent emission calculation cor- 
responding to two-gluon emission in the anti-fc^ algorithm which in the soft limit 
works like a perfect cone. At the two-gluon level we have four terms corresponding 
to the independent emission of soft gluons in the energy ordered regime xi ^ X2- 
These contributions are depicted in figure ^ The contribution to the squared matrix 
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element for ordered two-gluon emission is the same for each of the diagrams in fig- 
ure H, up to a sign. The double real (labelled (a)) and double virtual contributions 
(labelled (d)) can be expressed as 



{pi.ki){pi-k2){P2-ki){P2-k2) 



which in terms of the energy fractions Xi and X2 introduced in section |3.1| becomes 
simply 

M/(*a.) = 256^1^ 1^-^^^,^^^^-^. (5.2) 

Since the calculation that follows below is intended for highly collimated jets, i? <C 1, 
we shall take the small angle limit of the above result, 61,62 ^ 1. A similar result 
holds for the one-real one- virtual terms (b) and (c) in figure ^ with a relative minus 
sign. We are now in a position to compute the jet mass distribution at the two gluon 
level for the independent emission Cp term. 

We start by noting that the integration region for all graphs can be divided 
according to whether the real gluons ki and /c2 are inside or outside the triggered 
jet. We have four distinct regions: ki, ^2 both outside the triggered jet, ki, /c2 both 
inside the triggered jet or either of the gluons inside and the other outside the jet. 
The condition for a given gluon to end up inside or outside the triggered jet depends 
on the jet algorithm we choose to employ. In the anti-fcf algorithm the condition is 
particularly simple when considering only soft emissions; such an emission k is inside 
the jet if it is within an angle R of the hard parton initiating the jet, else it is outside. 



Given this fact let us consider how the various diagrams (a)-(d) in figure ^ 
combine in the different regions mentioned above. Since we are computing the jet- 
mass distribution dE/dp for a fixed jet-mass p, the pure virtual diagram (d) makes 
no contribution and hence we shall omit all reference to it in what follows. In the 
region where both emissions are in the jet we shall treat the sum of graphs (a)-(c). 
Where the harder emission ki is in the jet and ^2 is out, graphs (a) and (c) cancel 
since the real ^2 does not contribute to the jet mass exactly like the virtual k2. This 
leaves diagram (b) which gives zero since the in-jet gluon ki is virtual and hence does 
not generate a jet mass. Hence the region with ki in and k2 out gives no contribution. 

Now we consider k2 in and ki out. The contributions with k2 real (a) and (b) 
cancel as the graphs contribute in the same way to the jet mass. The diagram with 
/c2 virtual (c) cannot contribute to the jet mass as the real emission ki lies outside 
the jet. 

Hence we only need to treat the region with both gluons in and we shall show 
that this calculation correctly reproduces the result based on exponentiation of the 
single gluon result. The summed contribution of graphs (a) to (c) can be represented 
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Figure 4: Diagrams contributing to independent two-gluon emission from a hard parton 
line. 



as 



dp 



j d^W[5{p- xi9l - X2el) -5{p- xi9l) - 5 {p - X29l)] 



(5.3) 



where we wrote the contribution to the jet mass from an emission with energy fraction 
X and angle 9 with respect to the hard parton as 2x (1 — cos6') ^ x9'^. 

To compute the leading double-logarithmic contribution and show that it cor- 
responds to the exponentiation of the order double-logarithmic term one can 



write 5 {p — Xi9f — X29l) 
approximation 



as (p — Xi9f — X292) and make the leading-logarithmic 



e (p - xi9l - X29l) ^ e (p - xi9l) e (p - X29l) , (5.4) 

which allows us to make the replacement 

5{p- xi9j - X29l) ^6{p- xi9j) e (p - X29l) +1^2. (5.5) 
Doing so and using the explicit forms of W and the phase space in the small 
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angle limit we get 



^^2 (<^s\'^ f d9ld9ld(j)dXidX2 f n2\ f n2 \ , i o1 

e - el) (7?^ - el) e (xi - ^2) e(i - x^). (5.6) 

Carrying out the integrals we straightforwardly obtain 

2 



dS2 
dp 



which is precisely the result obtained by expanding the exponentiated double-logarithmic 
one-gluon result to order and differentiating with respect to p. Thus the standard 
double-logarithmic result for the jet-mass distribution arises entirely from the region 
with both gluons in the jet. Contributions from soft emission arising from the other 
regions cancel in the sense that they produce no relevant logarithms. 

We shall now argue that for algorithms other than the anti-/ct relevant single- 
logarithmic contributions shall appear from the regions which cancelled in the argu- 
ment above, although of course the leading double-logarithms are still precisely the 
same as for the anti-Zct case. An analysis of such miscancelling contributions is there- 
fore necessary for a resummation aiming at next-to-leading logarithmic accuracy in 
the jet-mass. The logarithms we compute below correct the one-gluon exponentiated 
result for the jet-mass distribution at the single-logarithmic level starting from order 



Let us consider the situation in, for instance, the kt algorithm. When both ki 
and /c2 are within an angle R of the hard parton both soft gluons get combined 
into the hard jet and this region produces precisely the same result as the anti-/ct 
algorithm, corresponding to exponentiation of the one-gluon result. Moreover, when 
both ki and k2 are beyond an angle R with respect to the hard parton there is no 
contribution from either to the jet-mass. However, when ki is beyond an angle R 
and k2 is inside an angle R the situation changes from the anti-Zct case. This is 
because in the kt algorithm when the two soft partons are separated by less than R 
in angle they can be clustered together. The resulting soft jet has four-momentum 
/ci -|- /c2, when we use the four- momentum recombination scheme, and lies essentially 
along the harder gluon ki. Thus when ki is beyond an angle R it can pull /c2 out 
of the hard jet since the soft jet fci + ^2 which replaces k2 lies outside an angle R 
of the hard parton. This results in a massless jet and hence such a configuration 
gives no contribution to the jet-mass distribution. In precisely the same angular 
region the virtual /ci, real /c2 diagram (b) (obviously unaffected by clustering) does 
however give a contribution whereas in the anti-Zct case it had cancelled the double 
real contribution (a). The graph with ki real and k2 virtual gives no contribution 
as before. Thus a new uncancelled contribution arises for the kt (and indeed the 
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Cambridge-Aachen) algorithm from the region where the two real gluons ki and /c2 
are clustered, which can be given by computing the ki virtual k2 real graph in the 
same angular region. 

We now carry out this calculation explicitly. We consider the angles Of, 6^ and 
as the angles between ki and the hard parton, k2 and the hard parton and fci 
and k2 respectively. Applying the kt algorithm in inclusive mode means constructing 
the distances ufOl, uj\Q\ and i>j\d\2 along with the distances (from the "beam") 
bj\J^ ^bj\B? ^ where, for the e^e~ case we consider here, the energy w plays the role 
of the kt with respect to the beam in a hadron collider event. Now since 6^ > R^, 
Q\ < the only quantities that can be a candidate for the smallest distance are 
00202 and Ci;|6'f2- Thus the gluons are clustered and k2 is pulled out of the jet if 
012 < O2 < R. Otherwise /c2 is in the jet and cancels against virtual corrections. 

We can then write the contribution of graph (b) of Figure ^ in the clustering 
region 

dp v27r/ J 9( 02 2vr Xi X2 

e {el - R') e {ej - © {R' - el) . (s.s) 

Using the fact that in the small-angle approximation relevant to our study 

el2 = el + 9l-29i92Cos(t), (5.9) 
integrating over xi and X2 and using t = — one obtains 



dp ^ ^ \2tx) p J 91 t 2n 

e (t - 1) {9l - R^) e {Apt cos^ (j)-9l)e {R^/p - t) . (5.10) 

Carrying out the integral over 9\ results in 

d ^cluster _ Ar^2((^s\'^l f dtdcj) f 4ptcos^ 

e (t - 1) e (4ptcos2 - R^) e {R^/p-t) . (5.11) 



" ^.Cluster ^ — - / In liliriUl int 

dp 2 ^ V27r/ pJt2'K V / 



Now we need to carry out the t integral for which we note t > max ^1, 



4/9 COS^ ( 

.2 



In the region of large logarithms which we resum one has however that R^ p 
and hence 4^ ^^^2 ^ > 1- This condition is reversed only when p ~ i?^ a region not 
enhanced by large logarithms and hence beyond our accuracy. 
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It is then straightforward to carry out the t integral and doing so and extracting 
the leading singular behaviour in p produces the result 



d ^ i,,tcr ^ _4^2 /^^y 1 In i / ^ In^ (2 cos 0) 6 f cos - ^ 
dp \2tt/ p p J vr \ 2 



_0.728C^f^Viln-. (5.12) 
V27r/ p p 



This behaviour in the distribution translates into an next-to-leading logarithmic 
a? In^ - behaviour in the integrated cross-section, which is relevant for resummations 

S p O 1 

aiming at this accuracy. As we mentioned before the above finding of single logarith- 
mic corrections generated by clustering has also been reported before for the case of 
gaps between jets studies [Q. Note however that the logarithms found there had 
coefficients that depended on the jet radius as a power of the jet radius starting at 
the B? level. In the present case however the presence of coUinear singularities near 
the boundary of a jet of small radius R ensures that there is no power suppression in 
R and hence the logarithms generated are formally comparable to those we aim to 
control here and indeed those resummed in Refs. Likewise the clustering will 

also generate leading logarithms in the Eq variable, which are again unsuppressed 
by any powers of R and hence ought to be controlled. Lastly we point out that 
the logarithms generated by independent emission and subsequent kt clustering were 
actually resummed in Ref. E3] and that possibility also exists here. 



6. Conclusions 

We would like to conclude by emphasising the main points of our study. Given 
the current interest in the study of jet shapes and substructure for the purposes of 
discovering new physics at the LHC, it is worth examining the theoretical state of 
the art when it comes to looking at individual jet profiles in a multi-jet event. A 
step in this direction was taken for instance in Refs. [||, 0. In the present paper we 
have noted 

• Observables where one picks out for study one or more jets in multi-jet events 
are in principle non-global. The non-global logarithms will arise at next-to- 
leading or single-logarithmic accuracy in the jet-shape distributions. If one 
studies jet events with a fixed multiplicity by imposing a cut Eq on hadronic 
activity outside the high-p^ jets, there are non-global logarithms involving the 
ratio of the shape variable p and the energy fiow i?o, as was first anticipated 
in Moreover, there are also non-global logarithms in Eq/Q where Q is the 
hard scale of the process. These logarithms are leading as far as the distribution 
in Eq for a fixed p is concerned. 



- 20 - 



In the limit of narrow jets i? — )■ 0, one may naively expect the non-global 
contributions to the jet-shape distributions to vanish with R due to the ap- 
parently limited phase-space available for soft emission inside the jet. Here 
we have pointed out that the non-global logarithms do not vanish in the small 
cone approximation as mentioned for instance in Ref . . One finds instead at 
small R an effect that is independent of R and arises from the edge of the jet. 
However, in the limit of narrow well-separated jets R^ <C (1 —cos 6ij), where 6ij 
is the inter-jet separation, one can simplify the non-global contribution. In this 
limit, owing to QCD coherence and the nature of correlated multiple soft emis- 
sions, one can regard the non-global logarithms to arise independently from 
the boundary of each jet up to corrections that vanish as R^/ (1 — cos 9ij). For 
a measured jet one picks up logarithms in 2EqR^ / (Qp) while for each unmea- 
sured jet one has logarithms in 2Eq/Q. The resummation of these logarithms 
yields a factor for each jet j, which is the factor computed, in the large- Ai'c 
limit, for the hemisphere jet-mass in e^e~ annihilation in Ref. P3, again up 
to corrections vanishing as -R^/ (1 — cosOij). 



• The overall size of non-global logarithms depends on the precise values one 
chooses for Eq/Q, R and p. However, broadly speaking, we find the contribu- 
tion not to vary significantly with Eq and to yield corrections of order 15 — 20% 
in the peak region of the p distribution. Choosing Eq/Q of order p/R^ elimi- 
nates the non-global contributions from the measured jet but steeply enhances 
the contributions from the unmeasured jet and it is not an optimal choice for 
reducing the overall non-global contribution to this observable. 

• We emphasise that the above observations are valid for the anti-Zcj algorithm 
in which our ansatz for resummation of jet shapes in an arbitrarily complex 
event is to correct the one-gluon exponentiation with a product of independent 
non-global factors from each jet. We further have emphasised that switching 
to algorithms other than the anti-fct gives relevant next-to-leading logarithms 
in the shape distribution as well as leading logarithms in the Eq distribution, 
even within the independent emission approximation. Thus predictions for 
observables such as the one discussed in this paper, in those algorithms are 
prone to more uncertainty than our current study in the anti-fcj algorithm at 
least until such logarithms are also resummed. 

We would like to stress that the general observations in this paper are of applica- 
bility in a variety of other contexts. For instance, the issue of threshold resummation 
addresses limited energy flow outside hard jets, along the lines of the Eq distribution 
here. The consequent non-global logarithms and the issue of the jet algorithm have 
not been addressed to any extent in the existing literature. The same issues crop up 
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in the case of resummation in the central jet veto scale for the important study of 
Higgs production in association with two jets. 

We hope that an awareness of the nature and size of the non-global contributions, 
the simplification that occurs in the small R limit and our comments about the sit- 
uation in other jet algorithms will help to generate more accurate phenomeno logical 
studies for these important observables at the LHC In particular in future work we 
shall address in more detail the role of soft gluon effects and especially non-global 
logarithms on QCD predictions relevant to new physics searches at the LHC. As an 
existing example of such studies in the context of Higgs physics and the filtering 
analysis we can refer the reader to Ref. [^]. We shall aim to provide similar studies 
in the context of other shape variables in the near future. We also note that a study 
of resummed jet shapes and profiles would constitute an interesting test of QCD with 
early LHC data and we shall also generalise the present work with this aim in mind. 
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